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rf = myjm\{i = 2, 3) {my and ms denote the masses of the vector and B meson, respectively). In 
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the penguin annihilation and the hard-scattering emission contributions are essential to understand 
the polarization anomaly, such as in the B —)• fK* and Bg —)• (jxf decay modes. We also compare 
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I. INTRODUCTION 


Exclusive Bq (g = u, d, s) meson decays, especially Bq —)■ VV modes, where V stands 


or a 


a 


light vector meson (p, K*, u, 0 ), have aroused a great deal of interest for both theorists 
and in experiments 1^. In contrast to the scalar and pseudoscalar mesons, vector mesons 
can be produced in several polarization states. Thus, the fraction of a given polarization state 
is an interesting observable, apart from the decay widths. Phenomenology of the Bq —)■ VV 
decays offers rich opportunities for our understanding of the mechanism for hadronic weak 
decays and CP asymmetry and searching for the effect of new physics beyond the standard 
model. In general, the underlying dynamics for such decays is extremely complicated, but 
in the heavy quark limit (m;, —)■ cx)), it is greatly simplified due to the factorization of the 
hadronic matrix elements in terms of the decay constants and form factors. Based on this, a 
number of two-body hadronic B decays had been calculated in this so-called naive factoriza¬ 
tion approach [u|. However these calculations encounter three major difficulties: (i) for the 
so-called penguin-dominated, and also for the color-suppressed tree-dominated decays, the 
predicted branching ratios are systematically below the measurements, (ii) this approach can 
not account for the direct CP asymmetries measured in experiments, and (iii) the predic¬ 
tions of transverse polarization fraction in penguin-dominated charmless Bq —)■ VV decays 
are too small to explain the data, in which large such fractions are measured. All these in¬ 
dicate that this factorization approach needs improvements, for example by including some 


more perturbative QCD contributions 12|. In the current market, there are essentialh 
approaches to implement perturbative improvements: QCD factorization (QCDF) [13 
perturbative QCD approach(PQCD) 15|, and the soft-collinear effective theory (SCET) 


iree 

141 , 

16|. 


All these frameworks have been employed in the literature to quantitatively study the dy¬ 
namics of the Bq —)■ PP,VP,VV decays, having light pseudoscalar (P)and/or Vector (U) 
mesons in the final states. 

In the Bq —)■ VV decays, as the Bq meson is heavy, the vector mesons are energetic 
with By — ''^ b / 2 . As the spectator quark (u, d or s) in the Bq meson is soft, a hard gluon 
exchange is needed to kick it into an energetic one to form a fast moving light vector meson. 
The theoretical picture here is that a hard gluon from the spectator quark connects with the 
other quarks of the four-quark operators of the weak iuteraetiou Q . The u^er.y.ug theory 
is thus a six-quark effective theory, and can be perturbatively calculated [17J. In contrast 
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to the other two approaches (QCDF and SCET), the PQCD approach is based on the kx 


factorization formalism 


18 


20| . The basic idea here is to take into account the transverse 


momentum kx of the valence quarks in the hadrons, as a result of which the end-point 
singularity in the collinear factorization (employed in the QCDF approach) can be avoided. 
On the other hand, the transverse momentum dependence introduces an additional energy 
scale leading to double logarithms in QCD corrections. These terms could be resummed 
through the renormalization group approach, which results in the appearance of the Su- 
dakov form factor. This form factor effectively suppresses the end-point contribution of the 
distribution amplitude of the mesons in the small transverse momentum region, making the 
calculation in the PQCD approach reliable. It is worth mentioning that in this framework, 
the so-called annihilation diagrams are also perturbatively calculable without introducing 
additional parameters j^, 22|. The PQCD approach has been successfully used to study a 


number of pure anni 

nr 

in experiments [8|, 122 


lila 


ion type decays, and these predictions were conhrmed subsequently 


- 25|. Thus, in our view, this method is reliable in dealing with the pure 
annihilation-type and annihilation-dominated decays as well. 




Several years ago, H. Y. Cheng and C. K. Chua updated j4,l5| the previous predictions 
for Bq VV decays in the QCDF factorization approach by taking the transverse 
polarization contributions into account, and using the updated values of the parameters in 
the input wave functions and the form factors. In the PQCD framework, although many 
studies of the two-body i?q-decays are available j7-9|, a reappraisal is needed for the following 
reasons: (i) In the previous studies, the terms proportional to = myjm\ {i = 2, 3)” have 
been omitted in the amplitudes, especially in the denominator of the propagators of virtual 
quarks and gluons. As we point out later, these terms do bring the earlier PQCD predictions 
in better accord in terms of the measured observables in some problematic cases, such as 
the B (f)K* and Bg —)■ 00 decays, (ii) recent progress in the study of the distribution 
amplitudes of the vector meson, especially for the 0 meson, undertaken in the context of 
the QCD sum rules, may significantly impact on some of the calculations done earlier, and 
(iii) Experimental data for some of the Bq —)■ VV decays, such as the branching ratio 
and the polarization fractions of Bg 00, are now available. In addition, we work out 
a number of observables, such as 0||, 0x, A^p, A0|| and A0x for the hrst time in 

PQCD. Among others, we revisit the B ^ p (a;)0 decay modes, the direct CP asymmetry 
of which could help us distinguish the PQCD and competing approaches. A related issue 
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is the large fraction of the transverse polarization observed in some of these decays. In 
the PQCD framework, penguin-annihilation contribution is the key to understanding this 
phenomenon. Especially, the chirally enhanced (S-P)(S-|-P) penguin-annihilation gives rise 
to large transverse polarizations. Together with the hard spectator-scattering contributions, 
this could help solve the transverse polarization puzzle in the penguin-dominated Bq —)■ VV 
decays. 

This work is organized as follows. In Sec. Ull we outline the framework of the PQCD 
approach and specify the various input parameters, such as the wave functions and decay 
constants. Details of the perturbative calculations for the Bg —)■ VV decays are presented 
in in Sec. uni and the various input functions are given in the Appendix. Numerical re¬ 
sults of our calculations are presented in Sec. IIVI and compared in detail with the available 
experiments and earlier theoretical works. Finally, a short summary is given in Sec. |Vl 

II. FORMALISM AND WAVE FUNCTION 

Our goal is to calculate the transition matrix elements: 


M cx (VVineffiBq) , 


( 1 ) 


with the weak effective Hamiltonian 'He// written as 261 



Here, Vuh{x) and Vth{x) = d, s) are the CKM matrix elements, Cj(p) are the effective 
Wilson coefficient calculated at the scale /i, and the local four-quark operators Oj (j = 
1,..., 10) are defined and classified as follows: 

• Current-current (tree) operators. 


Oi — {baUi3)v-AiUi3Xa)v-A, O 2 — {baUa)v-A{Ul3^/3)v-A, 


( 3 ) 


• QCD penguin operators. 



O5 — {baXa)v_A 


V+A, Oq — 


{baXp)v_A y^X^'gq'Jv+A, (5) 


<?' 
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Electroweak penguin operators, 


O 


7 


O 9 


q' 

3 - 

-{baXa)v-A'^eq'{q0q'^)v-A, 

q' 


Os = -(&aX^)y_A ^eg/(g^g^)i/+A, 
q' 

3 - 

Oio = 2^baXi3)v-A'^eq'{q0q'^)v-A, 
q' 


( 6 ) 

(7) 


with the SU(3) color indices a and (3 and the active quarks q' = (m, d, s, c). The left-handed 
(right-handed) current V A are dehned as 7/^(1 ± 75 ). Following |27|, we introduce the 
following combinations a, of the Wilson coefficients: 


Ol — C2 + O1/3, 02 — Oi -|- O2/3, 

ai = C^ + ^^± 1 / 3 , * = 3 , 5 , 7 , 9 / 4 , 6 , 8 , 10 . ( 8 ) 


In the perturbative approach to hadronic Bq decays, several typical scales are encountered 
with large logarithms involving the ratios of these scales. They are resummed using the 
renormalization group (RG) techniques. Standard model specihes the Wilson coefficients at 
the electroweak scale mw, the W boson mass, and the RG equations enable us to evaluate 
the dynamical effects in scaling the Wilson coefficients in Eq. ([2]) from to mi,, the 
6 -quark mass. The physics between the scale mi, and the factorization scale A/i, taken 
typically as Ah — -^/mfeAgcD, can be calculated perturbatively and included in the so-called 
hard kernel in the PQGD approach. The soft dynamics below the factorization scale Ah 
is nonperturbative and is described by the hadronic wave functions of the mesons involved 
in the decays Bq -A- VV. Finally, based on the factorization ansatz, the decay amplitudes 
are described by the convolution of the Wilson coefficients C(t), the hard scattering kernel 
H{xi,bi,t) and the light-cone wave functions ^Mi,B{xj,bj) of the mesons [28|: 


A ~ 


J dxidx2dx3bidbib2db2b3db3 


xTr [C{t)^B{xi,bi)^M2{x2,b2)^M3{x3,b3)H{xi,bi,t)St{xi)e , 


(9) 


where Tr denotes the trace over Dirac and color indices, 6 j are the conjugate variables of 
the quark transverse momenta ktr, Xi are the longitudinal momentum fractions carried by 
the quarks, and t is the largest scale in the hard kernel H{xi,bi,t). The jet function St{xi) 
coming from the threshold resummation of the double logarithms In^ Xi smears the end-point 


singularities in Xi 


29| . The Sudakov form factor e from the resummation of the double 
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logarithms suppresses the soft dynamics effectively i.e. the long distance contributions in 


the large-fe region 


30 


m. 


In the PQCD approach, both the initial and the hnal state meson wave functions are 
important non-perturbative inputs. For Bg [q = u, d, s) meson, the light-cone matrix element 
could be decomposed as [3^, l33j 


d^ze^’^-^0\qf^{z)U0)\Bg{PB,)) = ^ 


{fBg + MBg)'y5 


.w - 


\/20 


s. 


( 10 ) 


13a 


where n = (1, 0, Ot) and v = (0,1, 0^) are the unit vectors of the light-cone coordinate 
system. Corresponding to the two Lorentz structures in the Bg meson distribution ampli¬ 
tudes, there are two wave functions 4>Bq{k) and (j)Bg{k), obeying the following normalization 
conditions: 


= = (11) 

where Jb^ is the decay constant of the Bg meson. Due to the numerical suppression, the 
contribution of (pB is often neglected. Finally, for convenience, the wave function of B meson 
can be expressed as: 




( 12 ) 


with the light-cone distribution amplitude 

(pBgix, h) = NBgX^{l - x^) exp 






^ 2i 2 


( 13 ) 


where Nb^ is a normalization factor and ujg is a shape parameter. For B^{B^) meson, we 


use (jjq = 0.4 ± 0.04 GeV, which is c 


well known decay modes 


3,Q, 


etermined by the calculation of form factor and other 


34l |. Taking into account the small SU(3) breaking and 


the fact that the s quark is heavier than the u or d quark, we use the shape parameter 
ujg = 0.5 ± 0.05 GeV for the Bg meson, indicating that the s quark momentum fraction is 
larger than that of the u or d quark in the B^ or 5° meson j^. 

The light vector meson is treated as a light-light quark-antiquark system with the mo¬ 


mentum = My, and its polarization vectors e include one longitudina’ 
ei and two transverse polarization vectors e^, which are dehned in j?. 


_polarization vector 


35|. Up to twist-3, 
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the vector meson wave functions are given by 36[ |: 


= -^ [Mv^L(j>v{x) + + Mv(j)v{x)] 

= -^ [Mv^T(j)vi^) + ixf^^v^x) + Mv'ie^^p<^757^eyn^u‘^(/)y(a:)] , (14) 

for the longitudinal polarization and the transverse polarization, respectively. Here e^^vpa is 
Levi-Civita tensor with the convention = 1. 

The twist-2 distribution amplitudes are given by 


<tiv(x) = ^^x(l-x) 1+a[^.c7(()+a|,,c7(() 
<t>v(x) = ^^(1 -x) I + ayc7(i) + (i^vCl'‘^{t) 


( 15 ) 

( 16 ) 


iT') 

with t = 2a: — 1, and fy are the decay constants of the vector meson, which for V = 
p,u,K*,(j) are shown numerically in Table [Tl For the Gegenbauer moments, we use the 
following values [36|, l37j |: 


= 0, = 0.03 ± 0.02 (0.04 ± 0.03) , 

= 44^ = 0-15 ± 0.07 (0.14 ± 0.06) 44^ = 0 (0.20 ± 0.07) , 

=0.11 ±0.09 (0.10 ±0.08) . (17) 


For the twist-3 distribution amplitudes, for simplicity, we adopt the asymptotic forms 


4(x) = Mx) = 

= 1^(1 + *^). '^ i ((*) = ^(“‘)' 


( 18 ) 


TABLE I: Input values of the decay constants of the light vector mesons, taken from 


vector 

/v(MeV) 

/^(MeV) 

P 

216 ± 3 

165 ± 9 

UJ 

187 ± 5 

151 ± 9 

K* 

220 ± 5 

185 ± 10 

</> 

215 ± 5 

186 ± 9 
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FIG. 1: Leading order Feynman diagrams contributing to the —)• VV decays in PQCD 


III. PERTURBATIVE CALCULATION 

At leading order, there are eight types of Feynman diagrams contribnting to the Bg —)■ VV 
decays, which are presented in Fig{T] The first row shows the emission-type diagrams, with 
the first two contribnting to the nsnal form factor; the last two are the so-called hard- 
scattering emission diagrams. In fact, the first two diagrams are the only contribntions 
calcnlated in the naive factorization approach. The second row shows the annihilation-type 
diagrams, with the first two factorizable and the last two nonfactorizable. 

In the following, we shall give the general factorization amplitndes for these Bg —)■ VV 
decays. We nse the symbol LL to describe the amplitnde of the {V—A){y—A) operators, LR 
denotes the amplitnde of the {V—A)(y+A) operators and SB denotes that of {S—P){S+P) 
operators resnlting from the Fierz transformation of the {V — A) {V + A) operators. For 
the Bg VV decays, both the longitndinal polarization and the transverse polarization 
contribnte. The amplitndes can be decomposed as follows: 

-^(e2, es) = ■ 4*)^^ + {(^i^ua0n^v'^e2*°'e^*^)A^, (19) 

where is the longitndinally polarized decay amplitnde, A'^ and A^ are the transversely 
polarized contribntions, and is the transverse polarization vector of the vector meson. 

The longitndinal polarization amplitndes for the factorizable emission diagrams in Fig. (a) 
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and (b) are as follows: 


j-l rl/A 

jLL{LR),L ^ -87iCFMBfv2 / dxidxs / bidbib^db^CpBixi, h) {[{-1 + X3)4>3{X3) 

Jo Jo 

+r3{2x3 - l)(03(a:3) + 03 ( 0 : 3 ))] E^f{ta)Kf{xi,X3{l - rl),bi,b3) 

+2r3(j)l{x3)E^f{tb)Kf{x3,Xi{l - rl),b3,bi)} , (20) 


where r* = and Cp = 4/3 is a color factor. The functions he/, ta,b, and E^f can be found 
in Appendix El There is no (S' — P){S + P) type amplitude, as a vector meson can not be 
produced through this type of operators. In the PQCD approach, the traditional emission 
contribution is also the dominant one. Unknown higher order perturbative QCD corrections 
will influence the emission contributions as well as those from other topologies. At present, 
although the next-to-leading order (NLO) contributions have not been completed, the vertex 
correction has been done and is used to improve the predictions for the decays B 


and the i? —)■ vr form factors 


, which allows us to estimate the stability of the emission 


diagram in NLO. The results quoted below are based on the leading order calculations, but 
we also estimate the uncertainties from the partial NLO contributions based on the available 
results, as explained in Sec. HVl 

The last two diagrams in the first row in FigHJare the hard-scattering emission diagrams, 
whose contributions are given below: 

rl /l/A 




d[x\ 


bidbib2db2(j)B{xi, bi)4>2{x2) 

X { [{ X 2 - 1)03(2:3) + 02:3(03(2^3) - 03(2:3))] Eenf { tc ) henf { a , 01 , & 1 , ^2) 

+ [(2:2 + 2:3)03(2:3) -02:3(03(2:3) + 0*3(0))] Eenf{td)Knf{,a, I32MM)] , (21) 


J^enf = Iby-TrOFoM^ J d[x] I bidbib 2 db 2 (pBixi,bi) 

X { [o((2:2 - O - 1 )( 02 ( 2 : 2 ) 03 (o) - 02(2:2)03(0)) 

+(o + o - 1)(0*2(0)03(0) - 02(0)0*3(0))) 

+ {X 2 - 1 ) 03 (X 3 ) (02(2:2) + 02(0))] Eenfitc)henf{a, / 3 i,bi,b 2 ) 

+ [o((o - o)(0*2(0)03(0) + 02(0)03(0)) 

+ (o + 2 : 3 )( 02 ( 2 : 2 ) 03 (o) + 02(2:2)0*3(0))) 

+O 03 (O)( 02 (O) -0*2(0))] ^en/(O)hen/(a, 02 ,&l,& 2 )} , ( 22 ) 
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M 


SP,L 

enf 


d[x] 


•l/A 


bidbib 2 db 24 > B { xi , bi )( j ) 2 { x 2 ) 


X {[ 03 (x 3 )(a ;2 -X3 - 1 ) +r3X3(03(x3) +0*3(X3))] Ee„/( 4 )/len/(a, A, & 1 , ^2) 

+ [4>3{x3)x 2 + r3X3{4)l{x3) - 4)1 {x 3))] E^nf{td)Knf{a, P2,bl,b2)} . ( 23 ) 


The functions tc(d)-, Eenf, h^nf, a, Pi for the nonfactorizable emission diagrams are also listed 
in Appendix As is well known, the hard-scattering emission diagrams with a light meson 
(pseudoscalar or vector) are suppressed. This can be seen from the hgures (c) and (d), which 
are symmetrical. But, compared with the hgure (d), the anti-quark propagator in hgure (c) 
has an additional negative sign. As a result, the two contributions cancel each other. 

Figures (e) and (f) are the factorizable annihilation diagrams, whose factorizable contri¬ 
butions are listed below: 

/■i 

^LL{LR),L ^ SCfTt/bM^ / dX 2 dX 3 / 62 ^ 6263^63 

Jo Jo 

X { [(I)2{x2)h{x3){x3 - 1) + ‘Jr2r3(j)2{x2){x3(j)l{x3) “ (X 3 - 2)(j)l{x3))] 

■Eaf{Qhaf{ai,P,b2,b3) 

- [-X24>2{x2)4>3{x 3) + 2r2r3(t)l{x3){{x2 - 1)4>\{X2) + {X2 + l)cl)^2i^2))] 

'Eaf(tf)haf{a2, P,b3,b2)} , ( 24 ) 


r-1 r-l/K 

= -16Ci7’/B7rM| / dx2dx3 / 62 ^ 6263^63 

Jo Jo 

X { [ 2 r 203 ( 2 ^ 3 ) 02 (2^2) + r 3 {x 3 - l) 02 (a: 2 )(^3(^3) + 03(^3))] 

■Eaf{te)haf{ai, P,b2,b3) 

+ [2r302(a:2)03(x3) -F r2X2(j)3{x3){(Jl{x2) - 02 ( 3 ^ 2 ))] 

■Eaf{tf)haf{a2,P,b3,b2)} , ( 25 ) 

and the related scales and the hard functions listed in Appendix 

The last two diagrams in Fig.l are the nonfactorizable annihilation diagrams. The ex- 
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pressions for the corresponding amplitudes are as follows: 


■^anf = 16 y-CF 7 rM| j d[x\ J bidbib2db2(pBixi,bi) 

X { [r2r3{4>l{x2){(pl{x3){l -X2+ X3) + (j)l{x3){x2 + X3 - 1 )) 

+02(3^2) (03(^3) (1 -X2- X3) + (/> 3 (x 3 )(x 2 -X3 + 3 ))) 

-X 2 ^ 2 (x 2 )^ 3 (x 3 )] Eanfitg)hanf{,Oi, / 3 i, bi, 62) 

- [r 2 r 3 ( 02 (x 2 )( 03 (a; 3 )(l + X2 - X3) + 03 (a; 3 )(x 2 + X3 - 1)) 

+02(^2) (03(^3) (1 -X2- X3) + 03(0:3) (0:3 - 0:2 - 1))) 

+ (0:3 - 1)02(0:2)03(0:3)] Eanf(th)hanf{a, / 32 , bi, 62)} , ( 26 ) 

fo /■! 

■^anf = 16 y-C fvtMI j d[x] J bidbib2db2(j)B{xi, bi) 

X { [ 03 (o: 3 )( 02 (a; 2 ) + 02 (a; 2 ))r' 2 (o :2 - 2 ) 

- 02 (o: 2 )( 03 (o: 3 ) - 03 (o: 3 ))r 3 (o :3 + 1)] E an J (t^)(o^, /3i, 6i, 62) 

+ [ 02 (o: 2 )( 03 (a:^ 3 ) - 4 >lix 3 ))r 3 {x 3 - 1) 

- 03 (0:3) (02(2^2) + 02(a;2))?^2O:3] -Ea„/( 4 )^an/(«, 02, ^1, ^2) } , ( 27 ) 

fo /"i/A 

= IQy -CfttMb j d[x] J bidbib2db2(j)B{xi,bi) 

X { [r2r3(0*(o:2)(03(o:3)(l - 0:2 + 0:3) - 03(o:3)(o:2 + X3 - 1 )) 

+02(^2)(01(0:3)(0:2 + 0:3 - 1 ) + 03(0:3)(X2 - 0:3 + 3 ))) 

+ (0:3 - 1)02(0:2)03(0:3)] Eanf(tg)hanf{a, 0 i, bi, 62) 

- [r2r3(02(o:2)(03(o:3)(l + 0:2 - 0:3) + 03(o:3)(l - 0:2 - X3)) 

+02(0:2) (03(0:3) (X2 + X3 - 1) + 03(0:3) (X3 - 0:2 - 1))) 

-0:202(0:2)03(0:3)] Eanf{th)hanf{a, 02 , ^ 1 , ^2)} • ( 28 ) 

For the i?(s) —)■ VV decays, the transverse polarization amplitudes of the two factorizable 
emission diagrams yield: 

pI 

= 87rCFMyv2'r'2 dxidxs 61^6163^6305(0:1, fei) { [03 (0:3) 

Jo Jo 

ErsUxs + 2)03(0:3) - X303(x3))] Eef{ta)Kf{xi, ^3(1 - r^), 61, 63) 

+r 3 ( 03 (x 3 ) + 03 (o: 3 ))F;e/(f 6 ) 6 .e/(o: 3 ,Xi(l - r 2 ), 63 , 6 i)} , ( 29 ) 
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j-l 

jLL{LR),T ^ -STrC^M^/y^ra / dxidxs / &irf6i63d630B(Xi, &i) { [^3 (xg) 

Jo Jo 

+r3iix3 + 2)03 (x 3) - X303 (x 3))] Eef{ta)Kf{Xi,X3{l - rj), 6i, 63) 

+r3((/)3(x3) + 03(x3))Eey(4)/le/(x3,Xi(l - r^), 63,^1)} • (30) 

The transverse polarization amplitudes of the two hard-scattering emission diagrams 
fig.(c) and (d) are given below: 

fo ri/A 

J^enf^ = 16y-7rCi;^M|r2 / d[x] I bidbib 2 db 2 (pBixi,bi) 

X { [(1 - X2)03 (x3)(02(x 2) + 02(a^2))] £'en/(4)hen/(a, A, &1, ^2) 

- [2r3(x2 + X3) (02 (2^2)03 (^3) + 02(3^2)03(3^3)) 

-X203 (x3)(02(x 2) -|- 02(3^2))] -^en/(^d)^en/(«, 02, , ^2) } , (31) 

fo /•1/A 

J^enf = 16y J d[x] J 61^6162^620^(3:1, ^l) 

X { [(X2 - 1)03 (x3)(02(x2) + 02(312))] T/gn/(4)/^en/(«, 01, 61, 62) 

+ [2r3(x2 + X3)(02 (x2)03(x3) - 1- 02(X2)03(X3)) 

-X203 (x 3)(02 (x 2) + 02(312))] T/gn/(/^d)6.en/(«, 02, 61, 62) } , (32) 

[2 /•! ri/A 

J^enf = ISy-TrCi^M^ J d[x] J 6id6i62d620B(xi,6i)0^(x2) 

X { [r3X3(03(x3) - 03(X3)) - (pj{X3){rl{x2 - 1) - Xgr^)] 

‘-S'en/(^c)^en//3l; ^1; ^2) 

+ [^3^3(03(^3) - 03(^3)) + (Pl{x3){rlx2 + ^3X3))] 

-^en/(^c/)^en/(^5/^27 ^2) } 5 (2^) 

[2 /•! /"l/A 

J^enf = 16y-7rCFM| J d[x] J 6i(i6i62ci620ij(xi,6i)0^(x2) 

X { [r3X3(03(x3) - 03(X3)) - 03 (3I3) (^s (X2 - 1) + Xgr^)] 

f{tc)h 

en /(a, 01, 61, 62) 

+ [l’33:3(03(3l3) - 03(313)) + (Pl{x3){rlx2 - TgXg))] 

■^en/(/"d)6'en/(ci, 02, 61, 62)} , (34) 
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■^Inf = ISy- 7rC'FM^r2 J d[x] I bidbib 2 db 24 >B{xi,bi) 

X {[ 2 r’ 3 (a ;3 + 1 - X 2 ){(pl{x 2 )(pl{x 3 ) - 02(2:2)03(^3)) 

+ { X 2 - 1)03 (x 3 )( 02 (x 2 ) “ 02 ( 2 : 2 ))] ^en/(^c)^en/(«, 01 , , ^2) 

+ [a;20^(a:3)(02(x2) -02(2:2))] Eenf{td)henf{(X, (32MM)] , ( 35 ) 

■^Inf = ISy- 7 rCFM^r 2 / ci[a:] J bidbib 2 db 2 (j)B{xi,bi) 

X {[2r3(x2 - X3 - 1)(02 (x2)03(x3) - 02 (^2)03 (2:3)) 

+ {X 2 - 1)03 (a:3)(02(x2) -02(2:2))] Eenf{tc)henf{a, fdl,bi,b 2 ) 

+ [0:203(2:3) (02 (2:2) -02(2:2))] Eenf{td)henf{<y, ( 32 ,bi,b 2 )} ■ ( 36 ) 

The transverse polarization amplitudes for the factorizable annihilation diagrams are: 

/■i 

j^LLiLR),N ^ 8 CF 7 rfBr 2 r 3 M^ / dX 2 dX 3 / 62^6253^53 

X {[(^2(2:2)03(2:3) + 02 ( 2 : 2 ) 03 ( 2 : 3 ))( 2:3 - 2 )- 
2:3(02(2:2)03(2:3) + 02(2:2)03(2:3))] Eaf { te ) haf { ai , l 3 , b 2 M ) 

+ [(2:2 - 1 ) (02(2:2)03(2:3) + 02(2:2)03(2:3)) + 

(0:2 + 1 )(02(2:2)03(2:3) + 02(2:2)03(2:2))] Eaf{tf)haf{a 2 ,l 3 ,b 3 M)] , ( 37 ) 

r-l pl/A 

= 8 CF'irfBr 2 r 3 MB J dx 2 dx 3 J 62^6263^63 
X {[(02(2:2)03(2:3) + 02 (2:2) 03 (2:3)) (2:3 - 2 )- 
2:3(02(2:2)03(2:3) + 02(2:2)03(2:3))] Eaf{te)haf{ai, 0 , 63, 63) 

+ [(2:2 - 1 )( 02 ( 2 : 2 ) 03 ( 2 : 3 ) + 02(2:2)03(2:3)) + (2:2 + 1 ) 

•(02(2:2)03(2:3) + 02(2:2)03(2:2))] Eafitf)hafia 2 , 0 , 63, 62)} , ( 38 ) 

pi p 1 /A 

= IQCFTrfBM^ / dX 2 dX 3 / 62 ^ 6263(763 

Jo Jo 

X { [r203 (X3)(02(a:2) + 02(2:2))] ^a/( 4 ) 6 .a/(Q;i, 0 , 62, 63) 

- [r 302 (2:2)(03(2:3) - 03(2:3))] ^a/(7/)ha/(a2,0,63, 62)} . ( 39 ) 

From Eqs. fl 2 ^ and ( 1371 ) . we hnd that large cancellations between the two annihilation type 
diagrams ( (e) and (f)) take place, as a result of which they are highly power suppressed. 
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These two symmetric diagrams will cancel each other due to the relative negative sign 
introduced by the anti-quark propagator in diagram (e). This agrees with the naive argument 
that the annihilation contributions are negligible, especially for two identical final state 
mesons. For Eq. fl55]) . although the cancellations are not as severe as in Eq. fl2TD and Eq. (jSZD, 
the contribution is also highly suppressed being proportional to r 2 r 3 . From Eqs. fl2S]) and 
(l39|) . it is interesting to see that no cancellations or power suppression are involved. The 
chiral enhancement here is important to explain the large direct CP asymmetry, generated 
by the strong phase and the transverse polarization fraction in the penguin-dominated B 


decays 

SecHYl 


39l |. The chirally enhanced penguin annihilation contribution will be discussed 


m 


For the nonfactorizable annihilation diagrams ((g) and (h)), we get: 


M 


LL,N 

anf 


= ^Inf = 16a/ I d[x] 


» 1 /A 


bidbib2db2(j)Bixi, bi 


X {[- 2 r 2 r 3 {(j)^{x 2 )(j)l{x 3 ) -F 02 (^ 2 ) 03 (^ 3 )) 

- 0 ^(X 2)03 (X3)(r2(x2 “ 1 ) “ ^ 3 X 3 )] (tg)/?■„„/(«,/3i, 61 , 62 ) 

+ [(I)2{x2)(l)l{x3){rlx2 - rl{x3 - 1))] Eanf{th)hanf{a, f32,bi,b2)} , (40) 


anf 


-MfnT = 16 


» 1 /A 


J .in JO 


^ i d[x] J bidbib2db2(j)B{xi,bi) 

X {[ 2 r 2 r 3 (()) 2 (a; 2 ) 03 (x 3 ) + 02(2^2)03(^3)) 

- 0 ^(X 2)03 (X3)(r^(x2 - 1) +^ 3 X 3 )] Eanf{tg)hanf{a,( 3 i,bi,b 2 ) 

+ [(p2{x2)(p3{x3){rlx2+rl{x3 - 1))] Eanf{th)hanf{a, f32,bl,b2)} , 


(41) 


■^anf = --^anf = J d[x] J bidbib2db2(j)B{Xl, bi) 

X { [r 2(2 - X 2)03 (X 3 )( 02 (a: 2 ) + 02(2^2)) 

r3(2:3 -h 1 ) 0 ^(X 2 )( 03 (X 2 ) “ 03 ( 0 : 3 ))] F^an/(tg)han/(a, 01 , ^ 1 , 62 ) 

+ [^^3(2:3 - 1)0^(2:2)(03(2;3) - 03(2:3)) + 

r22:203 (2:3)(02(2:2) + 02(2:2))] Eanf{th)hanf{a,P2,bi,b2)} , ( 42 ) 

This completes the derivation of the various contributions in the i?(s) —)■ VV decays. We 
now turn to the presentation of our numerical results in the next section. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 


We start this section by listing the input parameters used in our numerical calculations. 
The vector meson decay constants have been summarized in Tabled! Other parameters, such 
as the CKM matrix elements, QCD scale (GeV), the masses (GeV) and the decay constant 
of the B(,^) mesons (GeV) and the corresponding lifetimes (in ps) are taken from the PDG 


review 


40| and are given below: 


Ag = 0.25 ± 0.05, Mb = 5.279, Mb, = 5.366, /b = 0.21 ± 0.02, fB, = 0.24 ± 0.03, 
Tb±/o = 1.641/1.519, tb, = 1.497, mfe(pole) = 4.8, 

Vud = 0.97427 ± 0.00015, K. = 0.22534 ± 0.00065, Kb = 0.0035llK;);[()(}t 
Kd = 0.008671°;E, K = 0.0404l°:°°Ji, K = 0.9991461°;°™, 

« = ( 8911^)1 7 = ( 6811 °)°. ( 43 ) 


With three polarization amplitudes, A^, , and A'^, the decay width is expressed as 


r(B,„ ^ VV) = JK [I |2 +2(1 p + I A'^ 7 ] . 

where the analytic formulas for the amplitudes A^, A^ and A'^ can be found in 
our convention, given in Eq. flT^ . the helicity amplitudes are dehned as follows: 


(44) 



. In 


Ao = Al = -A^, An = V2A^, Ax = V2A^, 


(45) 


B. 


where the dehnitions of Ao,||,x are the same as those in [8[. In this work, we also predict 
their relative phases 0|| = arg{A\\/A q) and 0x = arg{Aj^/A q). The polarization fractions 
/l,||,± dehned as 






Aq 

1 + 

+ 

CS| 

Ax 

2 


(46) 


In addition to the direct CP asymmetry parameters, we also evaluate the following observ¬ 
ables: 


= (/i -/l)/(/i +/i). Ahp={h-h)/(h + fi), A7| = (7| - 7)/2 . (47) 

Of these, A0x = (0± — is being worked out for the Bq —)■ VV decays for the hrst 

time. 
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TABLE II: Updated branching ratios, percentage of the longitudinal polarization f]^ and the trans¬ 
verse polarizations f±, relative phases, and the CP asymmetry parameters A^p and A^p in the 
B —>• Bg — Bg —)• (/)(/!) and Bg —)• K*^K*^ decays calculated in the PQCD approach. 


Modes 

Br{lQ-'=) 

fL{%) 

f± (%) 

(/'ll (rad) 

</>±(rad) 

B^> K*^4> 

q 8+4-9 


q-l q + lii.^ 


2 T4+U-23 

^•r+-o.i9 

Exp 

9.8 ±0.6 

48 ±3 

24 ±5 

2.40 ±0.13 

2.39 ±0.13 

B+ ^ K*+(l) 

ITSTfr” 

57.0+^:^ 



2 iq+9-^^ 

^•ra-0.20 

Exp 

10.0 ±2.0 

50 ±5 

20 ±5 

2.34 ±0.18 

2.58 ±0.17 

Bg (^(/> 


34.7tf;? 

31.6lli 

2 01 +0-23 
^•r'r_o,23 

Z.UU_Q 21 

Exp 

19 ±5 

34.8 ± 4.6 

36.5 ±4.4 ±2.7 

2.711[5|^ ± 0.22 


Bg ^ 

0-39lou? 

50.0t^;2 

24.2l|6 

i.yO_Q 22 

1 qK+0.21 

i.yO_Q 22 

Exp°' 

1.10 ±0.29 

51 ±15 ±7 

28± 11±2 

1.75 ±0.58 ±0.30 


Bg K*^K*^ 

r ^+3.0 
0-+_2,4 

qo q + 12.1 

30.0l|? 

9 1 9+0.21 
^•-'-"^-0.25 

9 1 p;+0-22 

Exp 

28.1 ±4.6 ±5.6 31± 12±4 

38± 11±4 




Afi^p{%) 

A^p(%) 

TTTO 

A(/)|| (rad) 

A(j)±{rad) 

B^ 

0.0 

0.0 

0.0 

0.0 

0.0 

Exp 


4±6 

-11± 12 

0.11 ±0.22 

0.08 ±0.22 

B+ ^ k*+4> 

—1 n+*^-is 
-‘-•'“'-0.26 

—0 60'’"®'^^ 

'“'•‘^-0.11 

—0 05"*"*^'^^ 
'“'•'“'^-0.33 

-0.01 

Exp 

-1±8 

17±11±2 

22 ± 24 ± 8 

0.07 ±0.2 ±0.05 

0.19 ±0.20 ±0.07 

Bg (l)(f) 

0.0 

0.0 

0.0 

0.0 

0.0 

Bg 

0.0 

0.0 

0.0 

0.0 

0.0 

Bg K*^K*^ 

0.0 

0.0 

0.0 

0.0 

0.0 



“The experimental results are taken from 


m 


For the charmless Bg VV decays, it is naively expected that the helicity amplitudes 
Hi (with helicity i = 0, —, +) satisfy the hierarchy pattern 


Hn : H_ : H. = 1 : 


A, 


QCD _ /Aqc'Z)x2 


(48) 


mb mb 

In the naive factorization approach, longitudinal polarizations dominate the branching ratios 
of B decays. In sharp contrast to these expectations, large transverse polarization of order 
50% is observed in i? —)■ K*(p, B K*p and Bg —)■ (pcj) decays, which poses an interesting 
challenge for the theory. This shows that the scaling behavior shown in Eq. fl48p is violated. 
In order to interpret this large transverse polarization many mechanisms have been proposed, 


such as the penguin-induced annihilation contributions 


form-factor tuning and even onset of new physics ^ 


42| . hnal-state interactions 43| . 
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TABLE III: Updated branching ratios (in units of 10 “®) of i? —)• VV decays calculated in the PQCD 
approach. For comparison, we also give the updated theoretical predictions in the QCD factoriza¬ 
tion (QCDF) approach 4 ] and the previous predictions in the PQCD approach^]. Experimental 


data are from the Particle Data Group 


40 ( 


Decay Modes 

Class 

This work 

QCDF 

PQCD(former) 

Exp 


C 

ri Q'y+U.lU+U.Ub+U.UU 
'-'■^'-0.09-0.04-0.01 

n q+i'5+i.i 
U.a_o.4-0.2 

0.9 ± 0.1 ± 0.1 

0.73 ± 0.28 

p+p- 

T 

9fi n+10-1+1-4+1-5 

9 c: r + l. 5 + 1.1 

^^■^-2.6-1.5 

35 ± 5 ± 4 

24.2 ± 3.1 

BO pO^ 

E,P 

ri /| o-hO. 15-|-0.09-|-0.01 
^•^'-'-0.12-0.08-0.01 

ri riO-|-0.02-|-0.36 
^•'-'^-0.02-0.00 

1.9 ± 0.2 ± 0.2 

< 1.6 

B^ —)■ UJU! 

C,P 

n crn-ro.2i-ro.o9-ro.o5 
^•'-’^-0.18-0.07-0.05 

ri ^-h0.9-|-0.7 
'^■'-0.3-0.2 

1.2 ±0.2 ±0.2 

< 4.0 

BO j^*0p0 

P 

Q oH“l-3-hl.l-|-0.0 
'^■'^-l.l-0.9-0.1 

A ^-h0.6-|-3.5 
^■'^-0.5-3.5 

5.9 

q 4+1-7 

BO K*+p- 

P 

Q ^-|-3.1-|-2.2-|-0.6 
^■^-2.8-1.9-0.9 

0 q+l-l-l-4.8 
°-’^-l.0-5.5 

13 

< 12.0 

B^ K*^uj 

P 

A 'y+ 2 . 1 + 1 . 6 + 0.2 

^•'-1.5-1.3-0.3 

9 C-+0.4+2.5 
^•'-’-0.4-1.5 

9.6 

2.0 ± 0.5 

BO 

P 

Q ~h0.13-|-0.10-|-0.02 
'-'■'^'^-0.11-0.09-0.03 

n ^-h 0 . 1 -|- 0.2 

'-'■’^-0.1-0.3 

0.35 

0.8 ± 0.5 

BO k*+K*- 

E 

n 9-1-1-0.09-1-0.03-1-0.01 
'-'■^-’•-0.09-0.05-0.02 

n 1 + 0 . 0 + 0.1 
^■-‘--o.o-o.i 

0.11 

< 2.0 

B° pV 

P 

n ni O+0.007-I-0.001+0.001 
U.U 10 _Q QQg_Q QQ 2 _ 0.001 



< 0.33 

B^ 

P 

n nin+0-005-i-o.ooi+o.ooi 
U.U1U_Q QQ4_Q QQ2_0.001 



< 1.2 

B^ (/.(/. 

P 

n ni o~l“3.003-|-0.005-)-0.001 
U.U1Z_Q qq2_0,004-0.001 


0 . 01891 °;“®} 

< 0.2 

B+ ^ p+pO 

T 

iq cr+5.0-r0.4+l.l 
-‘-'J''-’-3.9-0.7-1.0 

90 

^'J-U-l.9-0.9 

17 ± 2±1 

24.0 ± 1.9 

B+ P+LO 

T 

19 1 + 4 . 5 -rO.l+O.l 
-‘-^■-'■-3.7-0.4-0.1 

1 0 q+3.2+1.7 

17 ± 2±1 

15.9 ± 2.1 

B+ p+K*^ 

P 

Q q+3.8+2.7+0.3 
^■^-3.3-2.4-0.5 

Q Q-hl.2-|-3.6 
^■"^-1.1-5.4 

17 

9.2 ± 1.5 

B+ p^K*+ 

P 

1 + 2 . 5 -rl.3+0.3 
^■-’--1.9-1.3-0.5 

C r-h0.6-|-1.3 
^■^-0.5-2.5 

9.0 

4.6 ± 1.1 

B+ uK*+ 

P 

A rj + 1.7+1.3+0.3 
^•'-'-1.3-0.9-0.3 

q ri-h0.4-|-2.5 
'^■^-0.3-1.5 

7.9 

< 7.4 

B+ ^ x*+K*^ 

P 

n c^H"0*23-|-0.13-|-0.02 
'-'■^'^-0.19-0.12-0.02 

n ^-h0.1-|-0.3 
'-'■'^-0.1-0.3 

0.40 

1.2 ± 0.5 

B+ p+(j) 

P 

n riocH“0-013H“0.003-1-0.002 
U.UZO_Q q;^2-0.004-0.002 



< 3.0 


As pointed out in the context of QCDF [3(], after taking into account the NLO effects, 
e.g., vertex-, penguin- and hard spectator-scattering contributions, the effective Wilson co¬ 
efficients become helicity dependent. Including these effects, for some penguin-dominant 
modes, the constructive (destructive) interference in the transverse (longitudinal) amplitudes 
of i? —?■ VV decays makes the total transverse contribution comparable to the longitudinal 
one, and the transverse polarization fraction may reach as high as 50%. 
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TABLE IV: Updated percentage of the longitudinal polarizations of i? —VV decays calculated 


in the PQCD approach compared with the updated theoretical predictions in the QCD factorization 


are from the Particle Data Group 


(QCDF) approach and the previous predictions in the PQCD approach [ 3 ]. Experimental data 


40|. 


Decay Modes 

'i’his work 

QCDF 

PQCD(former) 

Exp. 


n 1 9H"^-'J4-|-U.iD-|-U.UU 
'^•-'-'^-0.02-0.01-0.00 

pi QQ-|-U.Ui5-|-U.Ub 

'^•^^-0.04-0.37 

0.60 

0.75 ±0.14“ 

5° ^ p+p- 

n Qc+o.oi+o.oi+o.oo 

U-yO-0.01-0.01-0.00 

n Q9+0-01+0-01 
'^•^^-0.02-0.02 

0.94 

0.977 ±0.026 


pi p*7-|-0.04-|-0.03-|-0.06 
'^■'^'-0.06-0.04-0.06 

c-^-O.ll+O.SO 

^•^^-0.25-0.36 

0.87 


ojoj 

n cr+O.07+0.04+0.06 
^■^'J-0.10-0.02-0.04 

r, q^+0.01+0.04 
^•^^- 0.01-0.20 

0.82 


5° ^ K*^p^ 

n /^c^~l“0.03-|-0.03-|-0.00 
^•'^^-0.03-0.04-0.00 

pi oq-|-0.00-|-0.60 

'^•'^^-0.00-0.31 

0.74 

0.57 ±0.10 

5° ^ K*+p- 

pi p*o~h0.04-|-0.03-|-0.02 
^•'^^-0.03-0.03-0.02 

n c:q+0-02+0.45 
'^•^'^-0.03-0.32 

0.78 


K*^u! 

n /^c;H"0.05-|-0.024-0.00 
^•'^'^-0.05-0.02-0.00 

n c:q+0.07+0.43 
^•^°-0.10-0.14 

0.82 

0.69 ±0.13 

BO ^ x*Oj^*o 

pi r qH“0.07-|-0.02-|-0.02 
'^■^°-0.08-0.02-0.01 

n c:9+0.04+0.48 
'^•'^^-0.07-0.48 

0.78 

0.80 ±0.13 

BO ^ K*+K*- 

~ 1.0 

~ 1.0 

0.99 


5° ^ pV 

n Qc+o.oi+o.oi+o.oo 
^•^^-0.01-0.01-0.00 




B^ ^uj(t) 

pi Qyi-|-0.02-|-0.01-|-0.00 
^•^■^-0.02-0.02-0.00 




B^ (/.(/) 

pi Q7-|-0.01-|-0.01-|-0.00 
^•^'-0.01-0.01-0.00 


0.65 


-)■ p+p° 

n Qo+0.01+0.01+0.00 
^■^°-0.01-0.01-0.00 

0 QiH+0-01+0-02 
'^•^'^-0.01-0.02 

0.94 

0.95 ±0.016 

-5- p+OJ 

pi Q7-|-0.01-|-0.00-|-0.00 
'^■^'-0.01-0.00-0.00 

0 q«+ 0-01+0-02 

'^•^'^-0.01-0.03 

0.97 

0.90 ±0.06 

B+ ^ K*+p^ 

pi 'yr H“0.03-|-0.02-|-0.02 
^•'^-0.03-0.03-0.02 

pi /27-I-0.02-1-0.31 
'^•'^'-0.03-0.48 

0.85 

0.78 ±0.12 

B+ ^ K*^p+ 

pi 7pi-|-0.03-|-0.04-|-0.00 
^•''^-0.03-0.04-0.01 

pi /I oH-0.034-0.52 b 
^•^^-0.04-0.40 

0.82 

0.48 ± 0.08 

B+ ^ K*+oj 

n c/f+0.06+0.02+0.04 
^•'^'^-0.06-0.02-0.03 

pi ^74-0.034-0.32 
'^•'^'-0.04-0.39 

0.81 

0.41 ±0.19 

pH- _, I<^*+ I<^*0 n 7/1H-0.03-1-0.02-1-0.01 

n ^ I\ I\ U./^_QQ 4 _QQ 3 _QQ 2 

r, .c+0.02+0.55 
^•^^-0.04-0.38 

0.75 

0.75 ±0.25 

B+ p+(t> 

n Qc+o.oi+o.oi+o.oo 
'-'•^^-0.01-0.02-0.00 





“This is from BABAR data 


4d| . The Belle’s new measurement yields 0.21 ^q'22 ± 0.13 


^This mode is employed as an input for extracting the parameters for B —>■ K*p decays in ref. 4|. 


In order to interpret the observed large transverse polarization fraction in the penguin- 
dominated B —)■ VV decays, e.g., B —)■ K*</), B —)> both the PQCD and the QCDF 

frameworks rely on penguin annihilation. However, in QCDF, the penguin-annihilation 
amplitude involves a troublesome endpoint divergence, which is fudged by introducing non- 
perturbative parameters. Hence, in QCDF, one can £t the existing experimental data on 
the branching ratios, and the CP asymmetries by adjusting the annihilation parameters 
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TABLE V: Direct CP asymmetries (%) in the B —>• VV decays and comparison with the predictions 


from QCDF Experimental data are from the Particle Data Group 


the data is from the ref. 


40|. For 5° ^ A:*o(+)pO(-), 


5C| 


Decay Modes 

p+p- 
BO pO^ 

^ Lju; 

BO j^*0p0 

BO ^ K*+p- 
BO K*^uj 
BO 

BO k*+K*- 

B^ pV 

B^ ^uj(j) 

5 ° ^ cjict) 

5+ ^ p+p° 

B+ -)> p+w 
B+ ^ K*+p^ 
B+ ^ K*^p+ 
B+ ^ K*+uj 
B+ ^ x*+K*^ 
B+ ^ K*+(j) 


I'his work QUDF Exp. 


7 n 7+2.y+U.8+3.8 

''^•‘-5.2-5.4-6.0 

on+iV+14 

'J'J-16-26 


p qo-|- 0.50-|-0.66-|-0.00 
'^•oo_o_59_o.3i_o.00 



rn ^+11.9+6.3+5.0 

q+2+51 


'yo y+6.7+2.6+3.3 
“''^•'-6.2-6.0-0.9 

qn+15+16 

0U_i4_i8 


Q Q+0.6+2.8+1.1 
^•^-0.6-2.8-1.0 

1 c;+4+16 
43-8-14 

-6±9±2 

ir\A r + 1.2+2.9 + 0.0 

^^•^-1.5-3.4-0.6 

qO+1+2 

^^-3-14 

21± 15±2 

r ^+0.3+1.2+0.8 
'^•'^-0.3-1.3-0.9 

00+9+5 

^3_5_i8 

45 ±25 

0.0 

— 14+4+6 
-'■+-1-2 


9Q O+2.0+6.4+4.6 
^^•°-5.7-9.5-4.7 

0 


0.0 



0.0 



0.0 



0.0 

n q+o+0.4 
U'O-O-0.5 


ri rir—0.03+0.05+0.00 
'^‘^'^-0.01-0.03-0.00 

0.06 

-5 ±5 

1 1 n+l.8+2.4+0.9 
-‘■■‘-•^-2.0-2.5-0.6 

0+1+3 

0-1-4 

-20 ±9 

cycy y+1.1+2.6 + 0.4 
'-1.5-2.5-1.2 

^q+6+12 

^0-2-28 

31± 13 

1 o+0.2+0.2+0.1 

■‘-•'^-0.3-0.0-0.2 

n q+0+2 
'^■'^-0-0 

-1 ± 16 

Q 1 +3.3+1.3+0.0 
^•-L-3.2-3.5-0.3 

16+^+'4 

30_4_43 

29 ±35 

no n+4.6+0.2+0.7 
^O.U_4 2_2.2-1.4 

lg+1+17 

-‘-3-3-34 


-1.0 

0.05 

-1±8 


B+ 


p+(t> 


0.0 


Pa and (pA, which reduces the predictive power of the theory. In contrast, in the PQCD 
approach, the annihilation type diagrams can be perturbatively calculated without intro¬ 
ducing any fudge factor (or parameter), which allows us to predict the direct CP asymmetry 
and transverse polarization. The large transverse polarization fraction can be interpreted 
on the basis of the chirally enhanced annihilation diagrams, especially the {S — P){S + P) 
penguin annihilation, introduced by the QCD penguin operator Oq A nice feature of 
the {S — P){S + P) penguin annihilation operator is that the light quarks in the hnal states 
are not produced through chiral currents. So, there is no suppression caused by the helicity 


19 











TABLE VI: Updated percentage of the transverse polarizations /_l(%), relative phases (/>||(rod), 
4>±{rad), Ai^ii (lO^^rad), A4>±{10~^rad) and the CP asymmetry parameters A^p{%) and A^p{%) 
in B —)• VV decays calculated in the PQCD approach. 


Decay Modes 

f± 

(p\\ 

</>+ 

^CP 

^CP 


A(pj_ 


45.9lg;2 

Y68™ 

2 81'*'0-04) 

88.9+720.7 

-11-6++ 

-98.9+^+ 

-105+?'^ 

p+p- 

2.421°;?^ 

q 1 9+0.06 
'^‘-'-^-0.06 

3 15+O.O6 

'3-40-0.05 

-2-05l°:i 

39.0I+ 

IO.2II? 

9.581+3 

BO pO^ 

16.7l|° 

q -iq + 0.17 

'J'-'-'J-o.ig 

q -iq+0.17 
'3-4<J-0.19 

26.6l}|« 

-6O.0I+? 

-87.8l?+ 

-98.4l?+ 

B^ uio 


q 9rj+0.25 

o.^U_0_20 

q 9-1 +0.24 
'^•^-'--0.22 

_3 75+11.8 

0-‘O_i6.2 

17-0+1? 

105l?+ 

1081??:? 

BO ^ K*^pO 

16.91?:" 

4 57+0-02 

^•'4 ‘-3.06 

+-OO-3.06 

3 54+4-20 
'3-0"1-1.07 

-7-7il?;i 

-0.12l?;?2 

0 22'’‘4-8® 

BO K*+p- 

15.6l|5 

q q-| +0.23 
0-0l_o.21 

q q(-» + 0.22 

O.OU_o_21 

23.8I5;? 

-50.91+ 

128+:? 

127+43 

4^‘-4.3 

BO K*^oj 

18.3l|6 

9 1 q+0.21 
"^•-'-^-0.20 

2 14+0-24 
^•4+-0.19 

1 45+4-44 
4-+0_i 62 

—8 52+®-*^^ 

9 90 + 1.79 

^‘^^-1.89 

-12-0+:? 


iq 7+4.0 
4^- ‘-3.6 

9 9^+0.20 
"^•^'^-0.16 

9 qi +0.19 

~ 0.0 

~ 0.0 

~ 0.0 

~ 0.0 

BO k*+K*- 

~ 0.0 

q q/|+0.08 
'^*'^^-0.06 

q q7+2.60 

'^•'^'-0.09 

0 02+^'^^ 
o-u^_o.oi 

7^ q+21-1 

56.4l?+ 

-125+^^® 

4^4’-2.0 

B° pV 

Z.OU_0 76 

q 7'(q+0.22 
'^*''^-0.31 

q 77+0.24 

'-0.27 

~ 0.0 

~ 0.0 

~ 0.0 

~ 0.0 

B^ 

9 7'o+1.08 
'^-0.86 

q 'y'y+0.20 
* '-0.28 

q 70+0.20 

'^•'^-0.25 

~ 0.0 

~ 0.0 

~ 0.0 

~ 0.0 

^ (pep 

U-U'3-0.02 

q 9(2+0.20 
'^‘^'^-0.14 

'3-OU_o.17 

~ 0.0 

~ 0.0 

~ 0.0 

~ 0.0 


n 45+0-08 

4’-+0-0.06 

Qg 

q 1 q+0.07 

3-4O-0.10 

0 002+®'^®^ 
U.UUZ_Q QQ3 

n q9+0.25 
'^•'^^-0.64 

-0-111+? 

—0 75+4‘-44 
44- ‘4’-0.45 

B+ p+u 

1 1 c+0-38 
-'■‘-'■^-0.29 

2 54+0-44 
^-'4+-0.15 

9 c;7+0-16 

9 ("19+0.69 
Z.UZ_o 74 

76.21??;? 

70.7l?+ 

83 5+47.3 
oo.a-19 7 

B+ ^ iy*+pO 

11 q+2.3 

1 54+4-44 
4-4’+-0.14 

1 54+4-43 
4-4’+-0.15 

11.3111 

-34.01+ 

-26.41??? 

97 q+158 
^/•O-4 0 

B+ ^ K*^p+ 

13.71?;^ 

1 SI +0.20 
-'-■^-*--0.18 

1 Q-| +0.19 
-*-•^^--0.18 

-0.36+1? 

0.98l+° 

4-4y-0.36 

— 1 34+4‘-41 
4-4’+-0.49 

B+ K*+oj 

17.2l|4 

9 1 q+0.23 
^■-LO_o 20 

9 1 q+0.22 
^•-*-^-0.20 

11-21+ 

-19.91+ 

-37.9+;? 

-38.7+:? 

B+ ^ x*+k*^ 

12-91^:1 

1 qfi+0.20 

4-4^°-0.17 

1 55+0.18 
4-4'44-0.19 

791 +2.54 
'•^-'■-2.50 

— 15 1+4.2 

44’-4-2.6 

20.2I+ 

28.4+:? 

.B+ p+(p 

2-36lJ:?« 

q 7(2+0.22 
'^*''^-0.31 

q 77+0.23 

'-0.27 

~ 0.0 

~ 0.0 

~ 0.0 

~ 0.0 


flip. As a result, the polarization fractions satisfy 

h ~ /|| ~ f±- (49) 

Thus, in the PQCD approach, the penguin-annihilation together with the hard-scattering 
emission diagrams can explain the large transverse polarization fraction measured in exper¬ 
iments. 

We present our numerical results for the branching ratios, direct CP asymmetries, and 
some other observables introduced earlier in the text, in Tables nxi The dominant topolo¬ 
gies contributing to these decays are also indicated in the tables through the symbols T (the 
color-allowed tree contributions), C (the color-suppressed tree contributions), P (penguin 
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TABLE VII: Updated branching ratios (in units of 10 of Bg —)• VV decays calculated in the 


PQCD approach. For comparison, we also cite the updated theoretical predictions in the QCDF 
approach ^ and the previous predictions in the PQCD approach 8|. Experimental data are from 


the Particle Data Group 


^1 


Decay Modes 

Class 

This work 

QCDF 

PQCD(former) 

Exp 

Bg — )• 

K*-p+ 

T 

c\A rj+iu.y+i.2+u.u 


op\ Q+8.2+1.4+i.2 

^'-'•^-6.2-1.4-l.l 


Bg^ 

^*0p0 

C 

ri /i^+0.19+0.11+0.00 
'-'■^'^-0.15-0.07-0.03 

I q+2.0+1.7 
-^"^-0.6-0.3 

n qq + 0.09 + 0.14 + 0.00 
^•O'J_0,07-0.09-0.01 

< 767 

Bg — > 

K*°uj 

C 

n qc;+0.16+0.09+0.04 
^"^^-0.14-0.08-0.08 

-I+1.5+1.3 
-^•■^-0.5-0.3 

n. q-i+0.10+0.12+0.04 
'-’•O-*--0.07-0.06-0.02 


Bg 

k*+k*- 

P 

c- 4+2.7+1.8+0.3 

y ^+1.0+2.3 
'•'^-1.0-1.8 

f. 7+1.5+3.4+0.5 
'J-'-1.2-1.4-0.2 


Bg 

p°(f) 

P 

pi Qo+0.15+0.03+0.01 
'-'■^'^-0.05-0.01-0.02 

r, -1 O+0.01+0.09 
^•-‘-0-0.01-0.04 

pi 00+0.09+0.03+0.00 
'-'•^'^-0.07-0.01-0.01 

< 617 

Bg 

ojcj) 

P 

pi 1'y+0.10+0.05+0.00 
'-'■-'-'-0.07-0.04-0.01 

n 1 c+O.44+0.47 
'-'•-'-‘^-0.12-0.04 

n -iR+o.oo+o. 10 + 0.01 

^•-‘-^-0.05-0.04-0.00 


Bg 

p+p- 

P 

1 r+0.7+0.2+0.0 
-'■•'-'-0.6-0.2-0.1 

n rc+O.04+0.73 
^•'JO-0.04-0.53 

1 pi+0.2+0.3+0.0 
-'-•'-'-0.2-0.2-0.0 


Bg^ 

p^p^ 

P 

pi +0.39+0.22+0.00 
'-'■ '^-0.24-0.14-0.00 

pi q^+0.02+0.36 
'-'•'^^-0.02-0.26 

(.,+0.12+0.17+0.01 

^•^-‘--0.11-0.10-0.01 

< 320 

Bg^ 

P^UJ 

E 

pi pipiQ+0.003+0.001+0.000 
U.UUy_Q QQ3_Q QQ2_0.001 

n nr)4+0.0+0.005 

U.UU + _Q Q_Q QQ3 

pi pipiy+0.002+0.001+0.000 
U.UU 1 -0.001-0.001-0.000 


Bg 

OJOJ 

P 

pi y| pv+0.16+0.10+0.00 
^•^'-'-0.18-0.10-0.01 

n 1 q+0-02+0.21 
'-'•-‘-^-0.02-0.15 

n qQ+0.09+0.13+0.01 
'-'•'3‘^-0.08-0.07-0.00 



TABLE VIIL Percentage of the longitudinal polarizations Jl in Bg —)• VV decays and comparison 
with the QCDF approach Q] and the previous predictions in the PQCD approach 8]. 


Decay Modes 

This work 

QCDF 

PQCD(former) 

Bg 

K*-p+ 

n Qc:+*i.ui+u.ui+u.uu 
'-'•^^-0.01-0.01-0.00 

n Q9+*i.ui+u.ui 

'-'•^^-0.02-0.03 

pi Qqy+U.UUi+U.ULhii+U.UUU 

W.yo ( -0.002-0.003-0.002 

Bg^ 

^*0p0 

pi ry+0.06+0.06+0.01 
^•^'-0.10-0.08-0.00 

n Qn+0-04+0-03 
^•‘^^-0.05-0.23 

n /|ct;+0.004+0.069+0.006 
^•^^^-0.003-0.043-0.009 

Bg 

K*°co 

n c:n+0.07+0.11+0.01 
^•^^-0.08-0.15-0.01 

pi Qpi+0.03+0.03 
'-'•^'-'-0.04-0.23 

pi r 09+0.003 + 0.035+0.023 

U.OOZ_Q QQ2_0.029-0.013 

Bg 

k*+k*- 

pi Q+0.13+0.03+0.05 
^•^^-0.09-0.03-0.06 

n c;9+0-03+0.20 
'-'•'-'^-0.05-0.21 

n /)qo+0.051+0.021+0.037 
U.^OO_Q Q4Q_Q q23_o.oi5 

Bg 

p^cj) 

n oR+0.01+0.01+0.00 
‘-’•°^-0.01-0.01-0.00 

pi oq+0.01+0.02 
'-'•°^-0.00-0.18 

pi oypi+0.002+0.009+0.009 
U.O' 11-0.002-0.003-0.004 

Bg ^ 

UJ(j) 

pi ^Q+0.08+0.08+0.02 

'-'•'^^-0.09-0.09-0.02 

n QK+o.oi+o.oo 

^•‘^'-’-0.02-0.42 

pi y^ y| q+0.000+0.054+0.009 
U.^^'3-0.075-0.061-0.004 

Bg^ 

p+p- 

~ 1.0 

~ 1.0 

~ 1.0 

Bg 

O 

O 

Q.. 

~ 1.0 

~ 1.0 

~ 1.0 

Bg 

p^oj 

~ 1.0 

~ 1.0 

- 1.0 

Bg 

OJOJ 

~ 1.0 

~ 1.0 

~ 1.0 


contributions), and E (LU-exchange annihilation contributions). Theoretical uncertainties 
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TABLE IX: Direct CP asymmetries (%) in the Bg —)• VV decays and comparison with the QCDF 


approach and the previous predictions in the PQCD approach 


Decay Modes 


This work 


QCDF PQCD(former) 


K*-p+ 

Q i+i.4+i.u+u.2 

^•-‘--1.5-1.2-0.3 

-'--'■-1-1 

0 9+1.U+1.2+0.4 

°'^-i.2-1.7-1.1 


7+6.4+10.5+7.5 

'^^•‘-5.9-16.0-7.9 

46+15+10 

+'1-17-25 

f.-, 0+3.2+17.1+4.4 

'jr-O-4 7-22.8-2.3 


'70 1 H“2.9-1-13.IH-S.l 
— ^o.i_2_2_7.4-8.3 

“^*^-15-6 

f.9 -1+4.8+19.7+5.5 
l’^-l-3.9-12.6-1.9 

k*+k*- 

O O+2.5+0.5+0.0 
°-°-8.9-2.9-0.2 

91 -1-1-1-2 
^J-_2-4 

Q Q-hO.4-1-3.34-0.3 
^*'^-0.7-3.6-0.2 

p^(p 

A o-|-0.6-1-0.6-|-l-2 

Oq + l.O + lO 
°’l-0.0-36 

-in 1+0.9+1.6+1.3 

l'l-l-0.9-1.8-0.5 

uj 4 > 

OQ rjH-l.3-1-0.5-1-3.4 
^o.U_3 2_2.3-5 .i 

q-1-3-1-20 

“^-1-15 

q n+0.6+2.4+0.6 
'i-'l-0.6-2.4-0.2 


0.0 

n 9-I-O.4-I-O.5 
'^•^-0.3-0.2 

0.0 

+ 

1 

9 Q+0.7+1.5+0.2 
^•^-l.l-1.3-0.2 

0 

9 14-0.2-1-1.7-1-0.1 

“^■-*■-0.1-1.3-0.1 

p^p^ 

9 Q+0.7+1.5+0.2 
^•^-l.l-1.3-0.2 

0 

9 -I +0.2+1.7+0.1 

^■l-0.1-1.3-0.1 

p'^uj 

n -I+1.0+1.9+1.2 
•‘--‘•■-'■-1.5-4.4-1.4 

0 

n n+O-7+2.7+1.0 
'l-'l-0.5-3.9-0.4 

coco 

q q+0.8+1.5+0.5 
'J-'J-1.0-1.4-0.2 

0 

9 n+0.1+1.7+0.1 
'^•'1-0.1-1.3-0.1 


Bs 

Bs 

Bs 

Bs 

Bs 

Bs 

Bs 

Bs 

Bs 

Bs 

Bs 


TABLE X: Updated percentage of the transverse polarizations f±(%), relative phases (;i)||(rad), 
4>^{va.d), A(/)||(10“^rad), A(;i>_i_(10“^rad) and the CP asymmetry parameters A^p and A^p in 
Bs —)• VV decays calculated in the PQCD approach. 


Decay Modes 

/+ 

<^11 


^CP 

/I-L 

+*-CP 

AcP\\ 

A4>± 

Bs ^ K*-p+ 

9 qi 

^•Oi- 0,21 

8 07+0-01" 

'i-'l‘-0.09 

CO 

0 

1 

p c 

b b 

00 a 

9 71 -hU.by 
' -*--0.72 

'i'l n+^'^'3 

o+-U_io.5 

12.4^1:? 

1 2 5+4-5 

l^'+-4.8 

Bs K*^p° 

22.51^4 

1 04+2-52 

i-i‘+-o.io 

1 00+2-53 

I'l’l’-O.lO 

1 7 ir-1-21.2 
-'-'•^-13.0 

99 ri-1-29.9 
^^■'^-31.4 

-31.51?"" 

-36.5l??2g 

Bs K*^io 

26.1+?;^ 

0 1 O-1-0.33 
^*-‘■^-0.28 

9 9q-l-0.32 
^•^O-o 27 

-5.99t23-52 


30.7l35;9 

36.51313 

Bs K*+K*- 

97- 74-5.2 

Z /./-70 

q F:q-l-0.33 
^•^^-0.25 

0 t;4+0.36 

’1-O+-0.24 

45.41.'!;; 

— 02 0+5-5 

03 7+11-1 

93.41}!^ 

Bs ^ pV 

8.89l?;80 

q 1 1 4-0.10 
'^•-‘-'‘--0.09 

q 9Q-1-0.09 
*^•^^-0.09 

q 97-I-I.O7 
-1.19 

-32.81^:^ 

—43 7+9-9 
+' 3 -‘- 9.5 

-63.91106® 

Bg — y ujcj) 


q qQ-1-0.20 

'^*'^^-0.17 

q qF:-l-0.30 
o.OO-o _23 

—2 24+5'5’^ 

^•^+-5.45 


— 36 7+12-5 

q9 74-I6.5 
'^"^•'-18.6 

Bs p~^p~ 

~ 0.0 

0 4n+0.04 
'l-+'l-0.04 

q 97-I-O.I6 
-0.15 

0.0 

30 '5+75-0 

9 07+0.44 
2'8t_o.59 

-27.4l|9 

Bs ^ 

~ 0.0 

0 4Q+0.04 
'i-+'l-0.04 

q 97-I-O.I6 
-0.15 

0.0 

30.5lJ|° 

9 Q7-I-O.44 

^■^^-0.59 

-27.4l|9 

Bs — >• p^co 

~ 0.0 

q /1O+0.04 
"i-^O-O.OS 

z.oo-o 22 

0.0 

27.91^1 

-9.30l^;| 

-30.4111? 

Bs — )• ww 

~ 0.0 

0 4Q+0.04 

'i-+'i-0.04 

q 97-I-O.I6 
-0.11 

0.0 

30.8tJt3 

9 71 4-0.42 
' -'--0.52 

-26.7llf 
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quoted in the tables are estimated from three sources: the first error quoted is from the input 
hadronic parameters, such as the decay constants of the initial Bq and the final vector-mesons 
and the parameters in the distribution amplitudes of the initial and final states, which can 
be found in sec. mi and Eq. fH5D . The second error arises from the scale uncertainties, char¬ 
acterized by Aqcd = (0.25 ±0.05) GeV and the variations of the factorization scales t (from 
0.8t to 1.2t) detailed in Appendix The scale-dependent uncertainty can be reduced only 
if the next-to-leading order contributions are known. The last error is the combined uncer¬ 
tainty in the CKM matrix elements and the angles of the unitary triangle. In Tables UTIIVII 
and 1X1 we have combined these uncertainties by adding them in quadrature and show the 
resulting uncertainty, due to to the space limitations in the Tables. 

We now discuss these results. For the branching ratios, the most important theoretical 
uncertainty is the first error caused by the nonperturbative input parameters. In the PQCD 
approach, the wave functions are the primary important input parameters and they heavily 
influence the predictions of the branching ratios, as also discussed in [Ml- We have adopted 
the new updated wave functions. While, for the direct CP asymmetry parameters, the 
dominant uncertainty arises from the second error, which is caused by the unknown higher 
order QCD corrections. From the definition: 

_ BR(B ^ f ) - BR(B ^ f) 


A dir _ 

^CP — 


BR{B ^ f) + BR{B /) 


1 MB ^ /) 

2 _ 

MB ^ /) p 

1 MB ^ /) 

| 2 ± 

MB ^ /) P 


(50) 

it is apparent that the wave functions of the initial Bq meson and the final vector mesons are 
overall factors, hence they drop out in the ratio and do not provide significant contributions 
in the estimates of the direct CP asymmetries. Direct CP asymmetries are proportional 
to the strong phases originated from the hard part, and the NLO QCD corrections will 
influence the strong phases significantly. Not having these corrections at our disposal, we 
can only estimate them by varying the scales. The resulting theoretical uncertainty is larger 
than the one from the wave functions, and we assume that the variation of the scales is an 
adequate account of the NLO corrections at this stage. 

For comparison, the updated results of the QCDF approachjd, 5| and the earlier PQCD 
predictions |7H9| are also presented. We have updated the PQCD computations in this work 
and the main improvements are: (i) Use of the updated vector mesons distribution ampli¬ 
tudes with new estimates of the Gegenbauer moments and decay constants, and (ii) the treat- 
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ment of the terms in the decay amplitude proportional to the ratio rf = myJrn^Bii — 2, 3). 
Since wave functions are the most important inputs in PQCD, their improved knowledge 
is expected to yield improved estimates of the branching ratios, polarization fractions, and 
other observables. We recall that in the earlier PQCD computations, r^^-dependent terms 
in the denominator of the propagators of the virtual quarks and gluons were omitted. From 
Appendix A, we hnd that, although their contribution is formally power suppressed, it can 
numerically change the real and imaginary parts of amplitudes and enhance the transverse 
polarization component, especially for the penguin-dominant decays. To quantify this, we 
have listed the amplitudes, branching ratios and transverse polarization fractions of the 
penguin-dominant decays —)■ Bg 00 and the tree-dominant decay B~^ —)■ 

with and without the r^-terms in Tabl dXIl We note that for the two penguin-dominant de¬ 
cays, the impact of the r^^-dependent terms in the amplitudes of the annihilation part, as 
well as in the imaginary part of the emission diagrams, is numerically signihcant. Taking 
the factorizable annihilation diagrams as an example, in the range near xs —)■ 1 or a ;2 —)■ 0, 
the nonzero rf contributes a non-negligible imaginary part. So by keeping the rf-terms, 
the branching ratios are reduced, while the transverse polarization fractions rise. The two 
main improvements go in the right direction in explaining the observed branching ratios and 
the large transverse polarization fractions in 5 —)■ JF *0 and Bg cl)4> decays in the PQCD 
approach. For the tree-dominant decay B~^ —)■ p~^p^, however, the effect on the traditional 
emission diagrams produced by the r^-terms is tiny, as expected. This is further discussed in 
Appendix A. Thus, the improved PQCD treatment presented here yields better consistency 
with the data. 

In Table UTl we list the current experimental measurements in B^{B~^) —)■ and 

Bg ^ K 0, —)■ iC K*^ and Bg —)■ 00 decays and compare them with our theoretical 
results worked out in this paper. These decays are all penguin-dominated and are measured 
with large fraction of transverse polarization. For these decays, the naive factorization 
approach predicts too small branching ratio^ factor of 2 ~ 3 jj], due to the small 


contribution from the penguin operators. In j7-9,l4^, the authors have studied these Bq —)■ 


VV decays, but those predictions are not in good agreement with the currently available 
experimental data. The primary task is to bring up the branching ratios and explain the 
polarization anomaly in these decays. In our update, we explain the bulk of the data. 
However, we note that for Bg —?■ K*^K and Bg ^ K 0 modes, our calculated branching 
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TABLE XI: Amplitudes (10 ^), branching ratios (10 and the polarization fractions (%) with 

(and without) the r^-dependent terms in the —)• Bg (l)4> and B^ —)• p'^p^ decays. 


Modes 



Br 

/l 


emission -3.3-|-0.67i -0.66-1-0.06i 
annihilation 0.32-1.6i -0.43-1-0.84z 

0.64-0.05Z 

0.42-0.83Z 

9.8 

56 


emission -3.0-0.09z -0.71-0.012z 
annihilation -0.42-1.95i 0.05-|-1.28z 

0.69-1-0.03z 
-0.ll-l.38z 

15 

70 

Bs (pcpirf) 

Bg 

emission -2.8-|-0.37z -0.60-1-0.lOz 0.60-0.08z 16.7 34.7 

annihilation 0.68-1.2z -0.53-1-1.Oz 0.53-1.0z 

emission -2.6-0.02z -0.64-1-0.03z 0.63-0.005z 26.6 45 

annihilation -0.04-1.8z 0.18-1-1.8z -0.15-1.7z 

- 5 - p'^p^{rf) 

emission 3.0-1-5.9z 0.28-1-0.33z 

annihilation ~ 0 ~ 0 

0.27-0.29Z 

r-^ 0 

13.5 

98 

5+ - 5 - p+p° 

emission 2.8-|-5.8z 0.12-1-0.33z 

annihilation ~ 0 ~ 0 

-0.ll-0.29z 
~ 0 

13.3 

99 


ratios are ( 5 . 4 II 4 ) x 10 ® and (0.39^o;4g) X 10 ® respectively, which are much smaller than 
the data, though they are compatible with the QCDF predictions x 10“® and 

( 0 . 37 l°;“l°; 24 ) x lO'^ respectively. 

In Table IIIIl we have given our estimates of the B —)■ VV branching ratios for different 
topologies. For the penguin dominant decay modes (indicated by P in the tables), our 
updated predictions basically agree with the QCDF predictions, except for —)■ K*^u. 

Due to the constructive interference between the penguin emission contributions and the 
penguin annihilation contributions, our prediction for this decays is almost twice as large 
as that of QCDF, and it also comes out larger than the current experimental data. As the 
experimental error is still large, we wait for consolidated date from Belle-II experiment. For 
the color-suppressed decay B^ —)■ p^p^, the calculated branching fraction in this work is 
( 0 . 27 lQ;QglQ;Q 4 lQ;Q 5 ’) x 10 “®, while BaBar and Belle obtained ( 0.9 ± 0.32 ± 0 . 14 ) x 10 “® 
and ( 0.4 ± 0.41q3) x 10“® respectively, with the current world average being ( 0.73 ± 
0 . 28 ) X 10 “®. Our result, within errors, agrees with the Belle data. Judged from the isospin 
triangle, since the decay rate of B^ —)■ is so small, the rate for the decay B^ —)■ p'^p~ 
ought to be double that of B~^ P^P^■ In experiment, however, within errors, these 
two rates are equal to each other, which is puzzling. Thus, the experimental situation is 
still in a state of flux. In Table IIVI discussed in more detail later, we show that for the 
B^ —)■ p^p^ decay the longitudinal polarization fraction is as small as 12%. As is well 
known, for the color-suppressed decays, the longitudinal polarization contributions from 
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two hard-scattering emission diagrams largely cancel against each other. What’s worse, 
the remaining longitudinal polarization contributions are nearly canceled by those from 
the annihilation diagrams. On the other hand, the chiral enhanced annihilation diagrams 
and hard-scattering emission diagrams provide a large transverse polarizable contribution. 
In the end, the —)■ is almost totally dominated by the transverse polarization 

component. In Table HVl we adopt the BaBar data [4^, but note that Belle has provided a 
new measurement 0.21^0 22^0.13 which supports our theoretical calculations. Thus, it 


is important to have a rehned measurement of the branching fractions and the longitudinal 
polarization fractions for B ^ pp to draw dehnitive conclusion. It should be noted that if 
the next leading order corrections are included, the branching fraction of 5 —?■ p^p^ might be 


55| . The previous 


enlarged while its transverse polarization fraction f± will become smaller 
PQCD estimates for the B^ —)■ p^u decay rate exceeded the current experimental upper 
bound. In this work, this branching ratio is now lower than the upper experimental bound 
but is about a factor hve larger than the QCDF prediction due to a near cancelation of the 
color-suppressed tree amplitudes. In the framework of PQCD, although the color-suppressed 
tree amplitudes also almost mutually cancel, the decay can get signihcant contributions from 
the annihilation type diagrams so that the decay rate comes back up and is not as small as in 
the QCDF prediction. This, together with some other predictions, provides an experimental 
check on these two competing frameworks. 

In Table. CVl we have given the fraction of the longitudinal polarization component, fi 
for B —)■ VV decays, where we have compared them with the available data, and also with 
the previous PQCD and QCDF 4] approaches. Of these, the predictions for the decays 
B —)■ (f)p{oj) are worked out for the hrst time. For the B^ —)■ p^oj decay, we predict the 
longitudinal polarization fraction as small as 67%, which is due to a signihcant transverse 
polarization component, f±, from the penguin annihilation diagrams. The fi for this decay 
is in agreement with the QCDF prediction ^ but is signihcantly less than the previous 
PQCD prediction (87%). From Table. IIVI one also sees that for B^ —)■ K*^uj, our estimate 
for the longitudinal polarization fraction is in excellent agreement with the experimental 
data. We note that, for B^ —)■ uco, our predicted longitudinal fraction is 66%, while the 
QCDF approach yields 94% where the longitudinal contributions highly dominate the 
amplitude. In QCDF, as in B^ —)■ p^u, the penguin annihilation contributions are also tiny 
in B^ COLO. In PQCD, together with the hard scattering contributions, the considerable 
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penguin annihilation contributions yield a different result. For the —>■ 00, the previous 

PQCD prediction of the longitudinal polarization fraction is 65% j^, while our updated 
longitudinal polarization fraction is given by ~ 1, which is confirmed also in 56|. For the 
—)■ K*^p^ mode, our result is larger than the data, while in the QCDF framework, the 
central value is the same as the data, as this mode is used to extract the input parameters |4|. 

In this paragraph we shall discuss direct CP-asymmetries in the decays B —)■ VV shown 
in Table |V] and their current measurements. Though none of the current experimental 
measurements for the CP asymmetries shown in Table |V] is conclusive, they are in accord 
with our theoretical calculations. This, in turn, implies that the dominant strong phases 
in these channel estimated in our approach are in the right ball-park. From Table. |Vl one 
also notes that the CP asymmetries are large for the penguin dominant decays, but they 
are small for the color allowed tree-dominant decays and almost pure penguin-dominant 
processes, such as and K*<p. For P° —)■ p^u decays, our prediction is about 60%, 

while that of QCDF is only 3%. In PQCD, since the emission diagrams nearly cancel 
each other, the annihilation diagrams provide the dominant contributions. As direct CP 
asymmetry is proportional to the interference between the tree and penguin contributions, 
the sizable interference makes the CP asymmetry parameter large, reaching 60%. For 
the B^ —)■ p^p^/uu modes, the large penguin contributions from the chirally enhanced 
annihilation diagrams, which are at the same level as the tree contributions from the emission 
diagrams, make the the CP asymmetry parameter as large as 70%. On the other hand, for 
pure annihilation type decay P° since there are no contributions from tree 

operators, it is natural to expect that the direct CP asymmetry is practically zero. In 
summary, the entries in Tables IIIIl HVl and |V] show that for these B —)■ VV decays our 
updated predictions are in good agreement with experiment, and, broadly speaking, are also 
in agreement with the QCDF predictions jd]. 

In Table IVH we give the predictions for the perpendicular polarization fraction, f±, the 
relative phases, 0||(rad), 0_L(rad), A0||(1O“^ rad), A0j_(lO“^ rad), and the CP asymmetry 
parameters A^p and for the B —)■ VV decays for the first time in the PQCD framework. 
These remain to be confronted with the data. In fact, these variables are already experi¬ 
mentally measured in five channels: B^{B^) —)■ iF*(iF*+)0 and B^ K*^(p, Bg —)■ K*^K*^ 
and Bg —)■ 00, which are shown in Table m Our results are in good agreement with the 
data. 
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We now discuss the results for the Bg —)■ VV decays. Since the initial and the dual 
state distribution amplitudes (DAs) are the most important input parameters in the PQCD 
approach, our predictions for Bg —)■ VV decays in Tables rvTTlIVIIIllIXlIXl are almost the same 


as the predictions in js], as the DAs we adopted here are similar to those used 


m 


8|, except 


for the DAs of the 6 meson. For the Bg —)■ 


decay, the central values of the branching 


ratio and the longitudinal polarization fraction estimated in B are 35.3 x 10"® and 61.9%, 
respectively. It is apparent that neither the branching ratio nor the polarization fraction 
are in conformity with the experimental data, posted as (19 ± 5) x 10“® and (34.8 ±4.6)%, 
respectively. With the updated DAs of (j) meson, the current predictions of all the observables 
listed in Table mi agree better with the data. This can be confirmed by the similar updates 
in Bg —)■ TT+TT” and B^ —)■ K^K~ decays Ml]- Also, due to the terms proportional to the 
ratio rf = in the denominators, which we keep, their influence is expected to be 

more pronounced, as is larger than the other light vector-meson masses. For the rest of 
the decay modes, the numerical values of the polarization fraction are basically consistent 
with the former PQCD predictions js]. 

From Tables lyTTl and IVIIIl we note that for the color-suppressed decays Bg ^ K piu), 
the branching ratios in PQCD are smaller than in QCDF by a factor of 3 due to the near 
cancellation of the hard scattering contributions. On the other hand, chirally enhanced 
annihilation and the hard scattering diagrams enhance the transverse polarization contri¬ 
bution, making it comparable to the longitudinal polarization one. For Bg ujcj), the pure 
emission mode, the {S — P){S + P) densities in the hard scattering diagrams also contribute 
a sizable transverse polarization component. For Bg , due to the large transverse 

polarization contribution from chirally enhanced annihilation diagrams, the longitudinal po- 
larization fraction is as small as 40%, which is similar to Bg —)■ K*^K . We also emphasize 
the measurements of the modes Bg ^ K p{oj) and Bg —)> uj(j) to distinguish among the 
competing dynamical models in the interpretation of the polarization anomaly. 

Direct CP asymmetries of Bg —)■ VV decays are listed in Table HXl We note that they are 
small for the penguin-dominant processes, since the interference between tree and penguin 
contributions due to the former are too small, which is opposite to the tree-dominant process 
Bg —>■ K*~p'^, which also has a small direct CP asymmetry. For Bg —)■ QCDF predicts 

about 83% for direct CP asymmetry with large charming penguin contributions. But in 
the framework of PQCD, it is only —4.3% because this mode belongs to the pure emission- 














type processes. Hence, measurement of direct CP asymmetry in this mode will help us to 
distinguish the PQCD and the QCDF approaches. 

As is well known, SU(3) symmetry relates a number of —)■ VV and —)■ VV 

processes, such as Bg —)■ K*~p'^ and B^ —)■ p'^p~. In the PQCD approach, presented here, 
this relation is well satished: 


B{Bs ^ 


(24.0 


+ 10.9+1.2+0.0^ 
-8.7-1.4-2.4 ) 


~ B{B^ p+p-) 


(26.0 


+10.1+1.4+1.5^ 
-8.1-1.4-1.2 J) 


(51) 


in units of 10“®. On the other hand, SU(3)-breaking in the decay rates for B K*(j) and 
Bg —)■ 00 is signihcant, as can be seen in Table ITTl In the PQCD approach, the SU(3)-breaking 
effects are caused by the differences between the initial and hnal state wave functions, such 
as the shape parameter ub and ub^, as well as the decay constants of the B and Bg mesons, 
along with the Gegenbauer moments and the decay constants of the hnal vector mesons. 
They conspire to yield a cumulative 60% SU(3)-breaking effect. Other SU(3)-breaking effects 
he in between these two cases, as can be numerically calculated from the entries in various 
tables presented here. 

t/-spin symmetry, relating a number of i?(s) —?■ hih 2 {hi are light mesons) has been advo¬ 


cated in the literature 


57l |. For VV decays, it has been studied in and checked 


against the explicit QCDF estimates, and seems to hold well. Since we have calculated the 
B and Bg decays to VV in this work in the PQCD approach, we also check the f/-spin 
symmetry in some representative decays studied in 


Acp{Bg ^ K*-p+) = -Acp{B^ ^ Jr+p 


_^B{B^^K*+p-)T{Bg) 

B{Bg K*-p+) t{B) ’ 


Acp{Bg ^ K*^p^) = -Acp{B^ ^ 
Acp{Bg —)■ p'^p ) = —Acp{B^ K*^K 

Acp{Bs —?■ K*^K* ) = —Acp{B^ p'^p 


0 . r..o.MB^^K*^p^)T{Bg) 


B{Bg K*0p^) t{B) ’ 
B{B^ K*+K*-)T{Bg) 

’ B{Bg ^ p+p-) t{B) 

B{B^^p+p-) rjBg) 
B{Bg K*+K*-) t{B) ' 


(52) 


Using these U-spin relations as well as the branching ratios, the lifetimes of B and Bg 
mesons and the direct CP asymmetries in B decays, we can get the relevant direct CP 
asymmetries in Bg decays. This can be then compared with the explicit calculations in the 
PQCD approach to check whether the f/-spin symmetry works well or not. We show this 
comparison in Table IXII( where the entries in the last two columns have to be compared 
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with each other. We hnd that, within the calculational errors, the ?7-spin symmetry works 
well in the direct CP asymmetries in the PQCD approach as well. 

TABLE XII: The direct CP asymmetries (%) in Bg —)• VV decays via [/-spin relation together 
with the direct PQCD prediction. The branching ratios of B and Bg decays are in units of 10~®. 


modes 

Br 

Acp{%) 

modes 

Br 

Acp{%m 

Acp(PQCD) 

BO K*+p- 

8.4 

od r+i.^H-^.y+u.u 

Bg ^ K*-p+ 

24.0 

-8.4 

Q 1 +i.4+i.u+u.i; 
"^‘-'■-1.5-1.2-0.3 

BO K*^pO 

3.3 

o q+0.6+2.8+1. 
^‘^-0.6-2.8-1.0 

Bg K*^p^ 

0.40 

72.3 

7+6.4+10.5+7.5 

'^^•'-5.9-16.0-7.9 

BO k*+K*- 

■ 0.21 

90 O+2.0+6.4+4.6 
^^•°-5.7-9.5-4.7 

Bg p+p" 

1.5 

-4.1 

9 q+0.7+1.5+0.2 
■^•^-1.1-1.3-0.2 

^ p+p- 

n n oqH“0'50H-0.66H-0.00 

ZU.U V7.00_Q 5g_Q 3^_Q QQ 

Bg ^ K*+K*- 

5.4 

-3.9 

o O+2.5+0.5+0.0 
°-°-8.9-2.9-0.2 


V. SUMMARY 


In this paper, we have reexamined the branching ratios, polarization fractions, relative 
phases, and direct CP asymmetries in Bg —)■ VV {q = u, d, s) decays in the PQCD approach. 
Compared to the previous PQCD calculations, the updated longitudinal and transverse 
decay constants as well as the Gegenbauer moments in the vector mesons wave functions 
have been adopted, which allows us to reduce the parametric uncertainties in the branching 
ratios and other observables. What concerns the predictions of the polarization fractions 
and their relative phases, we have kept track of the terms proportional to the ratio rf = 
m\rjm\{i = 2, 3), which have been ignored in some earlier estimations. In addition, we have 
studied the decay modes B —)■ p{uj)(j) that have not been explored before. For the observables 
f±, 0|h ^0||) we have provided the first PQCD predictions. So, 

this work updates and goes beyond what is already known in this approach. 

Our numerical results are listed in the Tables in the preceding section. For the well- 
measured decay modes, such a.s B ^ K*(l) and Bg —)• 00, the updated PQCD predictions for 
all the experimental observables fare better than the previous predictions in this approach, 
improving comparison with experiments. In addition, in many B{Bg) VV decays, our re¬ 
sults agree with the updated QCDF predictions jd, 5|, as well as with the experimental data. 
Yet, in some other cases, our predictions and those in QCDF differ and we have discussed 
some of these decays, such as B^ —>■ involving the annihilation contributions. 
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For the tree dominated B ^ pp processes, our results respect the isospin triangle re¬ 
lations, while the experimental data, taken on the face value, shows significant isospin- 
violation. Our estimated decay rate and the polarization fraction in —)■ are in good 

agreement with the Belle measurement, but not so compared to the BABAR data. This 
calls for a rehned measurement oi B ^ pp decays in the future. 

From the entries in Tables HVl and rvTTTl we note that our updated longitudinal polarization 
fractions are in good agreement with the data and the predictions in the QCDF approach % 
0] in some topologies. But for the color-suppressed decay modes, B^ —)■ B^ —)■ p^oo, 

Bg ^ K and Bg ^ K oo, the longitudinal contributions dominate the decay amplitudes 
in the QCDF approach, while in this work, the transverse polarization contributions are 
comparable to the longitudinal polarization contributions, and are even dominant in the 
amplitude for B^ —)■ This provides the possibility of distinguishing between these two 

approaches. 


Table and HX] list predictions of the CP asymmetry parameters,which agree with 


experimental data, wherever available, and, generally, also with the QCDF predictions j^, 


le 


a 


7*0 


n 


in some topologies. For the color-suppressed decays, B^ —)■ p^p^, B^ —)■ ojoj, Bg ^ K p' 
and Bg ^ K u, both PQCD and QCDF predict large direct CP asymmetries. But for 
po central value of the QCDF prediction is only 3%, while the prediction 

of this work is about 60% due to the large annihilation contributions. For B^ K*~^u 
and Bg —)■ pcp decays, which are almost purely dominated by penguin contributions, we 
predict very small CP asymmetries, but QCDF predicts them to be of orders 0.56 and 0.83 
respectively due to the charming penguins, which needs to be confirmed by experiments. 

Our predictions for many 5° —)■ VV decays basically agree with the previous PQCD 
predictions 8j. But for a few penguin dominant decay modes, for example, Bg 00, Bg 
K 0 and Bg —)■ K*^K , the improvements are signihcant, especially in the polarization 
fractions. 
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Appendix A: Related Hard Functions 


In this appendix, we summarize the functions that appear in the analytic formulas in the 
Section imi The hrst two diagrams in Fig{T]are factorizable emission diagrams, whose hard 
scales ta(b) can be determined by 

ta = max{Y/x3(l - tDMb, l/h, I/ 63 }, (Al) 


tb = max{^Xi{l-r‘l)MB, l/h, I/63}. (A2) 

The function h^f consists of two parts: the jet function St{x) and the propagator of virtual 
quarks and gluons. 

hef{xi,X3,bi,b3) = Ko{l3bi)[9{bi-b3)Io{ab3)Ko{abi) 

+0(63 — bi)Io{abi)Ko{ab 3 )] St{x 3 ), (A3) 


with a = and (3 = -^/x^x^Mb- The jet function in the factorization formulas can be 

given as|58|: 

21+2T(3/2 + c) 


St{,x) = 


[x{l-x)W 


(A4) 


+ c) 

with c = 0.4. In the nonfactorizable contributions, due to the small numerical effect, we drop 


the jet func 


diagrams 591. 


ion in the nonfactorizable emission diagrams and nonfactorizable annihilation 


The evolution factors Ef.f{ta) and Ef,f{tb) in the matrix elements are given by 

Eef{t) = a,{t)exp[-SB{t) - S3{t)]. 


(A5) 


The Sudakov exponents are dehned as 


SB{t) = s{xi^,bi'^ + ^ 


(A6) 


Si{t) = s{xi^,b^ + s (^{1 - Xi)^,b^ ^ ^ h/ 


(A7) 
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where the s{Q,b) can be found in the Appendix A in the Ref. 18|. Xi is the momentum 
fraction of “quark” in vector meson, with i = 2,3. 

For the rest of diagrams, the related functions are summarized as follows: 

4= max{y^(1 - r^)x5Xi Mb,^\ [(^2 - 1)(1 - ^1) + Xi)][rl + X3(l - rf)] | Mb, 

1 / 61 , 1 / 62 }, (A 8 ) 


td = max 


{\l{^-rl)x3XiMB,^\ [X 2 {rl - 1 ) +a:i)]x 3 (l - r|) \Mb, 


1/61,1/62}. 


(A9) 


Eenfit) = as{t) exp[-SB{t) - S2{t) - 53 ( 6 )]! bi=b3- 


(AlO) 


henf{oi, A, ^1, A) — [6^(62 — 6i)/o(a6i)A’o(«62) + 6^(61 — 62)/o(Q'62)iFo(®A)] 




X< 2^0 

Ko (AMB 62 ), 


'/IMB 62 ), A^< 0 ; 

(3f > 0, 


(All) 


with i = 1,2 and 


a = ^(1 - rl)x2,XiMB, (A12) 

Pi = [{.X 2 - l){l - rl) + xi)][rl + x^iil - rl)], (A13) 

Pi = ^ 2(^3 - 1 )+a:i)]a: 3 (l-r^), (A14) 

The hard functions and the scales for factorizable annihilation diagrams Fig.(e) and (f) are 

te = max{Q;iMB,/3MB, 1/62,1/63}, 

6/ = max{Q;2MB,/SMB, 1/62,1/63}, (A15) 

Eaf{t) = as{t) ■ exp[-S'2(A - S^it)], (A16) 


haf{ai,P,b2,h) 


= iPMBb 2 ) [eib 2 - 63)//^'^ («.Mb62) Jo (a.MBbo) 

+A63 - 62)//^'^ (a^MBbs) Jo {aiMBb2)] • ^(2:3), 


with 


ai = 

0:2 = 

P = 


\l^-x^{l-rl) 

^J{l-rl)[rl + X2{l-rl)\, 

V/[(1 - rl){l - X3)][ri + X 2{1 - r|)]. 


(A17) 

(A18) 

(A19) 

(A 20 ) 
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For the nonfactorizable annihilation diagrams, the scales and the hard functions are 


^9 

= ma.x{aM 

B, \/\f3i\MB, 1/bi, l/b2}, 

(A21) 

th 

= max{ q/Mb, l/6i, 1/62}, 

(A22) 

Eanfit) 

= as{t) ■ exp[-SB{t) - S2{t) - Ssit)] \ 6^=63, 

(A 23 ) 


i,&2) - y 

9{bi - b2)H^^^ (aMBbi) Jo {aMBb 

2) 


+e{b2 - {aMBb 2 ) Jo («Mb6i) 




(^^MlMBbi) , A < 0, 

(A 24 ) 


1 

Ko , A > 0, 

a = 

\/(l -^ 3)(1 

- r2)[r|-b 0:2(1 -r|)]. 

(A 25 ) 

/Si = 

1 - - rl){l - X 2 ) - xi][rl + xsil - rl)], 

(A 26 ) 

/S2 = 

(1--2)(1- 

X3)[xi - 0:2(1 - rl) - rg]. 

(A 27 ) 
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